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Abstract: We study three-dimensional Yang-Mills-Higgs theories with and without a 
Chern-Simons interaction. We find that these theories admit a rich spectrum of vortex 
sohtons carrying both a topological charge and a global fiavour charge. We further 
derive a low-energy description of the vortex dynamics from a gauged linear sigma 
model on the vortex worldline. 



Contents 



1. Introduction 



1 



2. Dyonic Vortices in li = 2 + 1 



2 



3. The Dyonic Vortex Worldline Theory 



6 



4. Example: A Single U{2) Vortex 



10 



1. Introduction 

The discovery of vortices carrying genuine non-Abelian flux [1, 2] lias led to a series of 
studies aimed at developing an understanding of the dynamics of these objects. Much 
of this work has focused on the role of vortex strings in d = 3 + 1 dimensions, with 
the associated Higgs phases and the issue of dual confinement being of special interest. 
It has also been suggested that non-Abelian vortices might have relevance for hot or 
dense QCD. For reviews, see [3, 4, 5, 6]. 

In this paper, we consider non-Abelian vortices in d = 2 -|- 1 dimensional planar 
theories. We are motivated by the possibhty that these sohtons could be realised in 
condensed matter systems. Studying theories both with and without a Chern-Simons 
interaction [7] , we uncover a rich spectrum of vortices carrying both a topological charge 
and a global flavour charge. We flnd that the dynamics of these dyonic objects depends 
intricately on the balance between the charges and the Chern-Simons coefficient. 

The plan of this paper is as follows. In section 2 we introduce the model of interest, 
then find the Bogomolnyi bound on the energy of vortex configurations and show that 
the dyonic vortex mass is different in the cases with and without a Chern-Simons 
interaction. In section 3, we review how motion on the vortex moduli space is affected 
by massive matter and Chern-Simons interactions, and we rederive our vortex mass 
results from the perspective of a = + 1 gauged linear sigma model. Finally, in 
section 4, we look at the implications of our results in the case of a single vortex 
in a [/ (2) gauge theory. We find that there is a variety of possible BPS and non-BPS 
solutions for the moduli space motion, including circular orbits and more exotic looping 
trajectories. 
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2. Dyonic Vortices in = 2 + 1 

We will begin by introducing the theory we are interested in, and highlighting some 
related models. Then we will write down the energy functional for our theory and show 
that it is subject to a Bogomolnyi bound. This bound gives the mass of a vortex. We 

will sec that the mass depends both on the topological charge and on a Noether charge 
associated with the global flavour symmetry of the theory. The form of this dependence 
is different in the cases with and without a Chern-Simons interaction. 

We will work with a U (N) Yang-Mills-Chern-Simons theory in d — 2 + 1 dimensions, 
with a real adjoint scalar and Nf fundamental scalars Qi with real masses rrii, i — 
l,...,Nf. If the theory also has the appropriate fermions then it exhibits J\f — 2 
supersymmetry (i.e. 4 supercharges), which determines the bosonic Lagrangian: 



For N > 2, the Chern-Simons coefficient k must be an integer so that the partition 
function is invariant under large gauge transformations. For the Abelian = 1 theory, 
there is no such constraint. Generically, the masses rrii break the flavour symmetry of 
the model from SU{Nf) to U{l)^f-\ 

This model is part of a large family of theories admitting vortex solutions. We can 
obtain various related theories by taking limits of the above Lagrangian. 

• When K, = and rrii = 0, we have a Yang-Mills theory with massless fundamental 
scalar fields. This theory admits non-Abelian vortices. These were originally 
introduced in [1, 2], and have since been studied in some detail (for reviews, see 



• When K = but rrii 7^ 0, we have a Yang-Mills theory with massive fundamental 
scalar fields. The vortices of the o? = 3 -|- 1 version of this theory were discussed 



• When K, but rrii = 0, we have a Yang-Mills-Chern-Simons-Higgs theory hke 
those studied in [10, 11]. This reduces to the Maxwell-Chern-Simons-Higgs theory 
of [12, 13] in the = 1 case. 




(2.1) 



[3, 4, 5]). 



in [8, 9]. 
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• When K ^ 0, rrii = and — >■ oo, the Yang-Mills term vanishes, and we can 
integrate out 4> to get the Chern- Simons- Higgs theory with sixth order potential. 
The Abelian version of this theory was first introduced in [14, 15, 16] and its 
non-Abelian generalisation has been studied more recently in [17, 18]. 

The review [19] covers many of the Chern-Simons models mentioned above. Vortices in 
Yang-Mills-Chern-Simons theories with several Higgs fields and no fundamental matter 
fields have been previously studied in [20, 21, 22]. 

Our theory has vacua featuring different amounts of symmetry breaking. Wc will 
work with the Higgs phase, which admits topologically stable vortices. This phase 
exists only when Nf > N; for simplicity, we set Nf = N from here on. In the Higgs 
phase, the vacuum expectation values of the fields are: 



where a — 1, . . . ,N is the colour index. The flavour and gauge symmetries of the 
theory are each fully broken in this vacuum, but a subgroup of their product sur- 
vives: f/(A^)gaugc X f^(l)flivour ~^ ^(l)diag " ^otc that WC Can always shift the vacuum 
expectation values to set rrii = 0. 

By a novel application of the Bogomolnyi procedure [23, 24, 12], we will now find a 
lower bound on the energy of the field configurations of the theory. Requiring that this 
bound should be saturated will give us a set of first order differential equations whose 
solutions are vortices. 

We first write down the energy functional H of the theory. This is the Noether charge 
associated with a combination of a time translation and a gauge transformation. The 
gauge transformation is chosen so that the infinitesimal variation of any field under the 
combination is proportional to that field's covariant time derivative. 





(2.2) 



H = 



/ d^x iTr {El + B') + iTr ((I^o^)^ + (^^a^)^) 



t _ 1^ 

' 27r 



V 
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Note that we define Ea = Foa and B = Fu. We now follow the usual Bogomolnyi 
procedure, and complete the square. 



H 



-I 



|P±g,|2 + Tr {\B ± \e {q,q^ - ^Kcf^ - v''))' 



^Tr 



20 (-^i^B + i [(Po^i) qi - qi (2^0%)^] + ^T^o^E^ + ^ [^o0, </>]) 



± (TrB) ± z J] (Po^i) - iP^qi)^ q^ 



rrii 



(2.4) 



In this expression we have used the notation V± = Vi± iT>2. 



If K and nii are chosen to be zero then we may set = 0, so the third line in (2.4) 
disappears. However, if either k or rrii is nonzero and we have nontrivial fields gj, Aq,, 
then the equation of motion for (j) includes source terms, so we cannot set = 0. Thus 
for the third line in (2.4) to disappear, we must invoke Gauss' law: 

+ ^ [(^ogO - qi (^^ogO^j + + ^ [^00, 0] = 0. (2.5) 

Of the two terms in the last line of (2.4), the first integrates to ±27rA;f^, where k is 
the topological charge. The integral of the second term can be written as ± Qifrii, 
where we define 



Qi = i J d^x (qi (Voqi) - (X>o?i)^ ^i) ■ 



(2.6) 



The Qi are conserved Noether charges associated with the f/(l)^^^ flavour symmetry 
of the theory. If we take the trace of Gauss' law (2.5) and then integrate over space, 
we find that J2i Qi = 

All of the terms in the first two lines of (2.4) are non-negative squares, so we see 
that the energy H of our configuration is subject to a Bogomolnyi bound, which we 
interpret as the vortex mass. 



H > 



27rkv^ + ^Q, 



rrii 



vortex • 



(2.7) 



This result is remarkable because it means that sometimes increasing the value of 
l^iQiTTiil may decrease the vortex mass. This is an unusual property: typically, 
increasing the size of a charge on a soliton will increase the soliton's mass. 
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To saturate the bound (2.7), a configuration must obey Gauss' law (2.5) and solve 
the first order Bogomolnyi equations that arise from setting the squared terms to zero 
in (2.4). Given a configuration that satisfies these constraints, it is straightforward 
to show that the quantity + Ylii Qi'^^i is strictly positive if A; > and strictly 

negative if A; < 0. For the case A; = 1, we will show in the next section that the 
minimum possible value of Mvortex for given rrii is 2nv'^ + mpK, where rup is the most 
negative of the mj values if «; > or the most positive of the values if k < 0. 

In the special case k = 0, the Bogomolnyi process works slightly differently. The 
absence of the kB term in Gauss' law means that in (2.4) we can make one sign choice 
for the first line and the topological charge term, and a separate choice for the second 
and third lines and the fiavour charge term. As a result, we have a stricter Bogomolnyi 
bound: 



Once again, the bound is saturated if and only if Gauss' law and the Bogomolnyi 
equations hold. 

We call the vortices in this theory dyonic because Myortex has contributions from 
two types of charge, the topological charge k and the Noether charge Ylii Qif^i- The 
significance of the difference between the bounds (2.7) and (2.8) is that increasing 

Qi''T^i\ always increases the vortex mass in the k — case, but may either increase 
or decrease the vortex mass in the k ^ case. 

The dyonic vortices we have considered in d = 2 + 1 dimensions are related to several 
other kinds of dyonic solitons in different numbers of dimensions. These include dyonic 
instantons in d = 4 + 1 [25], the original dyonic monopoles in d = 3 + 1 [26, 27, 28], 
semilocal vortices in = 2 + 1 [30]^ and the dyonic domain walls known as Q-kinks 
in c? = 1 + 1 [33]. A wide variety of soliton states, including vortices and their dyonic 
extensions, were considered in [34, 35, 36, 37]. 

It is interesting to compare the dyonic vortex masses we have calculated with the 
masses of other types of dyons. For dyonic monopoles and Q-kinks, the mass can be 

^Semilocal vortices can arise in Abelian gauge theories with multiple flavours of matter, and are 
related to sigma model Q-lumps. The moduli space of semilocal vortices suffers from non-normalizable 
zero modes corresponding to the scale parameters of the vortices [31, 32]. Our non- Abelian model 
avoids this problem because it contains equal numbers of colours and flavours: it hence has no vortex 
scale parameters, so the vortex moduli space is free of non-normalizable zero modes. 




(2.8) 
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written schematically as M = ^/Q^~+~Ql, where Qa and Qb are the two types of charge 
that make these sohtons dyonic [29, 33]. We can therefore regard dyonic monopoles 
and Q-kinks as bound states, each containing two distinct objects with masses Qa and 
Qfj. For dyonic instantons and semilocal vortices, the mass is schematically closer to 
M — Qa + Qb [25, 30], so it seems these objects are only threshold bound states. 

The dyonic vortices in our study, meanwhile, can be regarded as true bound states of 
two different types of object. The first object is a vortex with a topological charge but 
no flavour charges; this has mass 27r\k\v'^. The second object has flavour charges but no 
topological charge, so it is made up of excitations of the squark fields qf. Excitations 
of the about the vacuum have masses •\/e^(g")^/2 + (mo — m,:)^. The dyonic vortex 
mass Mvortcx is smaller than the sum of the masses of the constituent objects, since 
the (qf) part of the squark excitation masses does not contribute to Myortex- Hence the 
dyonic vortex is a classically stable bound state. 

3. The Dyonic Vortex Worldline Theory 

The moduli space approximation is a useful tool for understanding the low-energy 
dynamics of vortices [38]. For the theory given by (2.1), with k and set to zero, slow 
motion on the vortex moduli space is described by a nonlinear sigma model. In this 
section, we will review how this motion is modified when we reintroduce nonzero rrii 
and K. We will then reformulate our description of the motion as a gauged linear sigma 
model, and see how from this perspective we can rederive the results of the preceding 
section by applying the Bogomolnyi procedure to the d — + 1 energy functional. 

In the theory (2.1) with = and k = 0, the set of static physical solutions of the 
Bogomolnyi equations with topological charge k forms a moduli space Ai. This has 
collective co-ordinates and Kahler metric gpq, where p,q = 1, . . . , 2A^|A;|. Low-energy 
vortex dynamics is described by geodesic motion on A4: 

L = i(7,,(X)XPX'^. (3.1) 

The metric has an SU{N) isomctry arising from the unbroken 5'f/(A^)diag remnant of 
the gauge and flavour symmetries of the theory. We can flnd a set of N Killing vectors 
/Cj associated with the isomctry such that Yli^i ~ ^■ 

If we allow rrii ^ but keep k, = 0, then a potential is induced on the moduli space. 
This potential can be written in terms of the Killing vectors ICi [8, 9, 39]: 

N 

Vra = Yl (^^^n K^i) ^ 9pr (3.2) 



6 



If we instead allow k 7^ but keep rrii = 0, then a magnetic field JF e fi^ {/A) is induced 
on the moduli space. Writing this locally as ^ = dA and working to leading order in 
Ke/v, we find the dynamics on A4 is modified by an extra term in the Lagrangian 
[10, 13]: 

L« = -KAp{X)XP. (3.3) 

If we allow both rrii and k 7^ 0, then the Lagrangian for motion on Ai includes 
both the Chern-Simons term (3.3) and a potential similar to (3.2). We might expect 
that the potential (3.2) is modified by terms of order rriiK,; we have not calculated this 
from the nonlinear sigma model, but in the next subsection we will check that it is 
true, at least for /c = 1, by using a gauged linear sigma model. 



The Gauged Linear Sigma Model For One Vortex 

For an alternative description of motion on the one- vortex moduli space, and to clarify 
what happens when k and rrii are both nonzero, we now assemble a standard d = + 1 
gauged linear sigma model, with a U{1) gauge field and J\f = (0,2) supersymmetry. 
This 'worldhne theory' for a non-Abehan vortex is straightforward to solve, but it is 
more interesting than the corresponding theory for an Abelian vortex because non- 
Abelian vortices have an internal orientation as well as a spatial position. If we ignore 
the spatial position of our non-Abelian vortex and set rrii — 0, then the moduli space 
of vortex configurations is M. = CP^~^, with homogeneous co-ordinates (pi G C, i = 
1, . . . ,N. These co-ordinates give N fundamental scalar fields in the worldline theory, 
and roughly speaking they correspond to the orientation of the vortex in the U{N) 
gauge group of the parent theory. The fields (pi are subject to the D-term constraint 

N 

£> = J]99i99/-r = 0, (3.4) 

i=l 

where r is a constant which is determined to be 27r/e^ [1, 2, 40]. In the case = 0, 
K — 0, the Lagrangian (3.1) is L = \'Dt(pif, where Vtipi — dt(pi — iAcpi and A is the 
gauge field for the U{1) worldline gauge symmetry. The scalars cpi enjoy an SU{N) 
fiavour symmetry corresponding to the SU{N) isometry of the metric on M.. 

If we introduce nonzero masses for the Qi of the parent theory, then the gauged 
linear sigma model is deformed by the introduction of masses rrii for the (pi [8] : 

L = \Vtipif -^(p,^a - TTiif cpi. (3.5) 
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Here, a is an adjoint scalar field in the M = (0, 2) gauge multiplet. The deformation 
(3.5) is consistent with the J\f = (0, 2) supersymmetry of the theory, but for generic 
rui it breaks the flavour symmetry from SU{N) to U{1)^~^. When we integrate out a 
and A in (3.5), we recover the moduli space potential (3.2). 

If we now introduce a Chern-Simons interaction into the parent theory, then a cor- 
responding term is induced in the gauged linear sigma model [10]: 

L = \Vtipif - ^ ifii {a - rriif Lpi- k{A + a) . (3.6) 

i 

Note that we can again integrate out a and A to find the potential on }A: 

y = - ^(m, - mjf\ipi\^\ipj\^ + - ^mi|</?j|^. (3.7) 

i<j i 

The first term is independent of k, and matches the potential in [8]. The second term 
is the promised order uiiK modification. Note that the first term is non-negative, so V 
can be minimised (subject to the constraint (3.4)) by choosing p e 1, . . . , such that 
rUpH < rriiK, then setting 

{r if i —p ^ ^ 

(3.8) 
otherwise. 

Prom the Lagrangian (3.6), we can compute Gauss' law, 

N 

i J] {(Vm) - (Av'i)^) = « (3.9) 

i=l 

^^Qi^i^ , where = i (^(^(^i) - (^(^i)^) . (3.10) 

i 

The Qi are conserved Noether charges corresponding to the f/(l)^~^ flavour symmetry 
of the worldline theory. As we will see soon, they are related to the Noether charges 
Qi (2.6) of the parent theory. 



The Dyonic Vortex Mass 

We are now ready to apply the Bogomolnyi procedure to rederive the results of section 
2. As in the d = 2 -|- 1 case, we can find the energy functional for the theory as the 
Noether charge associated with a combined time translation and gauge transformation. 

Hid = \T^m? + {o- - rriif + Ka. (3.11) 
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When we complete the square, we get 



Hid = \T^m - - mi)^i 



+Tr 



a \ n — I 



(3.12) 



Imposing Gauss' law (3.9) gives the lower bound Hid > Qiirii. Note that if rrii = 0, 
then the vortex has mass 27rv^. The energy Hid gives corrections to this mass which 
arise from the fields on the worldline, so we can write 



rrii 



(3.13) 



This formula is remarkable because it shows from the worldline perspective that appro- 
priately chosen values for the charges Qi can give a negative contribution to the mass 
of the vortex. This is an unusual property for a soliton, but it matches what we found 
in the parent theory in section 2. 

By eliminating the auxiliary fields a and A from (3.11), we can see that the minimal 
possible value of Mvortex for given rrii is achieved by setting (pi = and then choosing 
(pi to satisfy (3.8), so that V is minimised. We then have Myortcx = 27rf^ + mpK,, where 
rUp is the most negative of the rrii values if /t > or the most positive of the rrii values 
if K < 0. This minimal value of Mvortex is guaranteed to be strictly positive by the 
fact that the vacuum expectation value (gp given in (2.2) must be strictly positive to 
permit the winding of the vortex to lie in the pth flavour, as demanded by (3.8). 

In the special case k, = 0, the Bogomolnyi process works slightly differently, just as 
in three dimensions. In particular, we have a choice of signs when we complete the 
square: 



Hid = \1^m ± ^(cr - mi)(pi 



±Tr 



=FEa 



rrii. 



(3.14) 



This gives us a stricter bound on the energy once we apply Gauss' law. 



H.d> 



rrii 



Mvortex = \2-Kkv 



rrii 



(3.15) 
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Our results (3.13), (3.15) match the k = 1 versions of the corresponding results 
(2.7), (2.8) from the parent theory if we make the identification Qi = Q^. Note that 
under this identification, there is agreement between the constraints Qi — nk and 
X^j Qi — 1^, and the right hand side of (3.13) is necessarily positive. Just as we saw in 
section 2, increasing Qiiriil always increases the vortex mass in the k — case, but 
may either increase or decrease the vortex mass in the k ^ case. 

If we are interested in a multi-vortex system, it is again possible to provide a gauged 
linear sigma model for the dynamics, this time by using a brane construction [1]. 
However, the metric induced on the moduli space when we quotient the gauged linear 
sigma model by the action of the gauge group no longer matches the metric Qpg of the 
nonlinear sigma model. Consequently, the fc-vortex gauged linear sigma model does 
not give a strictly accurate description of fc-vortex moduli space dynamics, although it 
may still be useful for finding qualitative properties of the motion. 

4. Example: A Single U{2) Vortex 

In this section, we will apply our results to the case of a single vortex in a C/(2) Yang- 
Mills- Chern- Simons theory including two flavours of matter with masses m and — m. 
We will obtain the effective potential on the moduli space of vortices, and use it to 
study some features of the dynamics. In particular, wc will find that the gauged linear 
sigma model perspective makes it straightforward to find BPS and non-BPS solutions 
corresponding to several different types of moduli space motion. 

If m = 0, the moduli space is = CP^. When we switch on the mass m, a 
potential given by (3.7) is induced on M.. Using co-ordinates 9 e [0, tt] and cf) e [0, 27r) 
on S^, we have 



The first term describes a standard sigma model on S^. The second term is the Dirac 
monopole connection, as was explained in [10], and the final two terms give the potential 
induced by m. Changing the sign of m corresponds to changing which pole we identify 
as north. Changing the sign of k corresponds to performing a parity transformation 
9 ^ 71 — 9. In what follows, we choose both m and k, to be non- negative without loss 
of generality. 

The theory described by (4.1) has a conserved angular momentum: 



L=- 9'^ + sin^ 9 (f)'^ +- {cos 9 - 1) 4> - rm^ sin^ 9 - Km cos 9 . 



(4.1) 



J ^r(j) sin^ 9 + K,cos9 ^ Qi - Q2. 



(4.2) 
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There is also a conserved Hamiltonian which we can write as 



// = 1^2 ^ V^s(9), where V^s(0) = — — ^T^^ + rm'^ sin^ + mK cos ^. (4.3) 

4r sin ^ 

The shape of the effective potential V^ff determines the dynamics of the system. Writing 
c = cos 9, s — sin 9, we find 

Vi,m = 2™^.- (c^ - + ^ - 1) (c3 - . (_L ^ 1) + £_) . (4.4) 
The Shape Of Kif When k = And m^O 

When J — 0, Ves has a maximum at 9 — 7r/2 and two minima at 9 — 0,7r. When 

I J| > 0, Ves — > 00 as ^ — >■ or TT. For < | J| < 2mr, V^s has a maximum at ^ = 7r/2 
and two minima at the values of 6 given by sin^^ = ^| J|/(2mr), as shown in figure 1. 
For I J| > 2mr, V^g has a single minimum at 6* = 7r/2. 

The energy of a solution sitting at a minimum of Ves is given by 

i 2rm'^ \J/(2mr)\ for I J| < 2mr , , 

' rm^ [(J/(2mr))2 + 1] for | J| > 2mr. 



Veff(^) 




To understand this, we turn to the Bo- 
gomolnyi equations that arise from (3.14). 
These are satisfied by setting ^ = along 
with either 9 — 0,9 — ironj) — 2m sign(mJ). 
For |J| < 2mr, solutions to the Bogomol- 
nyi equations are precisely the configurations 
that sit at minima of V^g; these saturate the 
Bogomolnyi bound i/icj > QiiTLi\ = \mJ\. 
For \J\ > 2mr, the Bogomolnyi equations ^ 3^ ^ ^ 

cannot be solved, and the energy exceeds the 424 
Bogomolnyi bound even at the minimum of Figure 1: A sketch of the effective poten- 
VeS- tial when k = and < | J| < 2mr. 

The Shape Of When k>0 

For nonzero Ves can take on a rich variety of different shapes. We find different 
patterns of behaviour in the cases < k < mr/2, mr/2 < k < Amr and 4mr < k, and 
in this section we will consider each of these cases in turn. 
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In every case, we have Ves — >■ oo as 6' — or tt, except where stated otherwise. 
The Bogomolnyi equations arising from (3.12) are solved by setting ^ = and either 
^ = 0, ^ = 7ror0 = 2m; note that by (4.2), is a function of 9. To satisfy ^ = 0, a 
configuration must sit at an extremum of Ves. We will refer to a minimum of Ves as 
'BPS' if it is located at a value of 9 such that ^ = 0, ^ = 7ror0 = 2m. Otherwise, 
we will say the minimum is 'non-BPS'. At BPS minima, the Bogomolnyi bound is 
saturated, so Vcs = fnJ . At non-BPS minima, V^s > itiJ- We will see that minima of 
Ves often occur at 6* = /±, where we define 



f±{hi, J) = cos 



I 2 I 2mr y V2mr/ 



J 



2mr 




We first describe the shape of Ves in the 
case < K < mr/2 for different values of J. 
When I J| < K, Veff has a non-BPS minimum 
and a maximum at values of 9 less than tt/2, 
and a BPS minimum at 9 = f^. A sketch 
is given in figure 2. When k < J < 2mr + 
fi;^/(8mr), V^s has two BPS minima at = 
f±, with a maximum in between. A sketch 
of Ves would then look like figure 1, but with 
the positions of the extrema shifted sideways; 
an exception is the special case J = k, when 
one of the minima is actually at ^ = 0. When 
J > 2mr + K^/{8mr), Vcs has a single non- 
BPS minimum at a value of 9 in (0, 7t/2). 








3£ 

4 



Figure 2: A sketch of the effective poten- 
tial when < K < mr/2 and \ J\ < k. 



Now suppose we decrease J Tintil it is negative. When J — —n, Kff has a BPS 
minimum at ^ = tt, and a non-BPS minimum and a maximum at values of ^ in (0, 7r/2). 
For 2mr(3[fi;/(4mr)]^/^ — 1) < J < —k, Ves has a non-BPS minimum at a value of 9 
in (7r/2,7r), and a non-BPS minimum and a maximum at values of 9 in (0,7r/2). A 
sketch of Ves would look similar to figure 2. Finally, when J < 2mr(3[fi;/(4mr)]^/^ — 1), 
Ves has a single non-BPS minimum at a value of 9 in (7r/2,7r). This concludes our 
description of possible shapes of V^s when < k < mr/2. 

The pattern of possible shapes is slightly different if wc choose mr/2 < k < 4mr 
because now —k, < 2mr(3[ft;/(4mr)]^/^ — 1). For —k, < J < 2mr(3[fi;/(4mr)]^/^ — 1), 
Ves has a single BPS minimum at ^ = /_; note that f- — n ior J — —k. For 
J < —K or J > 2mr(3[fi;/(4mr)]^/^ — 1), the behaviour matches what we had in the 
case < K < mr/2 for J < 2mr(3[«;/(4mr)]^/^ — 1) or J > —k respectively. 
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In the case k > Amr, Ves always has a single minimum. For \J\ < k, the minimum 
is BPS and is located aX d — f^. In particular, the minimum is at ^ = when J — k 
and at ^ = TT when J = —k. For J < — the minimum is non-BPS and is at a 
value of 9 in (0, 7r/2); for J > k, the minimum is non-BPS and is at a value of 9 in 
(7r/2,7r). Note that the results of this paragraph still apply in the special case m — 
if we redefine /_ = cos~^(J/«;). 

Possible Types Of Motion 

Each minimum of Vgs gives a stable solution of 
the equations of motion of the system. Usually, 
such a solution is a circular orbit around with 
fixed nonzero at a fixed value of 9. However, if 
the minimum occurs at ^ = or tt, then it gives 
a static solution sitting at a pole of S^. Also, 
for m — and \J\ < k, the BPS minimum of 
Vcs corresponds to a static solution with = 
sitting at a value of 6 in (0,7r): at this value of 
9, the effects of J and k, effectively cancel one 
another out. 

If we give the system an energy greater than Figure 3: Possible moduli space paths 
that required to sit at a given minimum of Vcs, ^^^^ nonzero k and m. 
then 9 oscillates around that minimum and the 

motion on can take two distinct forms, as shown in figure 3. For \J\ > k, the sign of 

always matches that of J. We then have motion similar to the lower, wave-like path 
in figure 3. Things are more interesting for \J\ < k. Then 

sign = sign ^ — cos 

If the minimum of V^e that we are oscillating around is at the critical value 9 — 
cos~^{J/k) then we always have motion similar to the upper, looping path in figure 3. 
If the minimum of that we are oscillating around is not at the critical value of 9, 
then we can have a wave-like path or a looping path depending on whether there is 
enough energy for 9 to cross the critical value as it oscillates. 
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